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Abstract 

A Gaussian ansatz for the wave function of two-dimensional harmoni- 
cally trapped anisotropic Bose-Einstein condensates is shown to lead, 
via a variational procedure, to a coupled system of two second-order, 
nonlinear ordinary differential equations. This dynamical system is 
shown to be in the general class of Ermakov systems. Complete in- 
tegrability of the resulting Ermakov system is proven. Using the ex- 
act solution, collapse of the condensate is analyzed in detail. Time- 
dependence of the trapping potential is allowed. 

PACS numbers: 03.75.Fi, 05 45. -a, 32.80.Pj 



1 Introduction 

We consider the mean-field theory for two-dimensional interacting Bose- 
Einstein condensates, described by the Gross-Pitaevskii equation |]. We 
suppose a time-dependent harmonic trap. Introducing a Gaussian ansatz for 
the amplitude of the wave-function of the Bose-Einstein condensate, there 
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results a coupled nonlinear system of second-order ordinary differential equa- 
tions for the time evolution of the system Usually, this system is stud- 
ied only in the isotropic case for which the condensate has equal width in all 
directions 0. The purpose of this work is to show that in the case of two 
spatial dimensions the isotropy assumption is not needed, since the nonlinear 
system for the Bose-Einstein condensate can be cast in the form of a com- 
pletely integrable Lagrangian Ermakov system || |7|, ||. Moreover, using the 
formalism of Lagrangian Ermakov systems, we are able to extract detailed 
information about the time evolution of the condensate. These results have 
an increased importance in view of the recent experimental achievement of 
Bose-Einstein condensation in quasi-one and quasi- two dimensions P|-jl2 



For instance, the present theory provides an accurate estimate for the critical 
time of collapse of the condensate (see equation (|33|) bellow), when collapse 
can occurs. 

Ermakov systems Jl3| have attracted interest in the last three decades or 
so due to both their physical applications and nice mathematical properties 
. The central feature of Ermakov systems is the existence of a constant of 



motion, the so called Ermakov invariant 14 , which generalizes the Lewis 



invariant for the time-dependent harmonic oscillator [I5[| . The Ermakov in- 
variant allows the construction of a nonlinear superposition law [HJ between 
the solutions of the differential equations of the Ermakov system. Moreover, 



the existence of the Ermakov invariant allows to linearize [17, 131 the Er- 



makov system in many cases. Other theoretical studies on Ermakov systems 



concerns its Lie symmetry structure [IS, 20], the existence of additional con 



stants of motion and the extension of the Ermakov systems concept to 
higher dimensions ||19|| . From the physical point of view, Ermakov systems 
have found applications in several problems, such as cosmological particle cre- 
ation ||22|| , nonlinear optics and propagation of shallow water waves 
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As will be seen in the following, two-dimensional Bose-Einstein conden- 
sates can be described by a nonlinear dynamical system of Ermakov type, ad- 
mitting a Lagrangian formulation. From the Ermakov property, one constant 
of motion, the Ermakov invariant, is immediately derived. The existence of 
a variational principle suggests the use of symmetry tools for the search of 
additional invariants. Here, we look for Noether point symmetries, which pre- 
serves the action functional up to addition of a numerical constant. We show 
that the Lagrangian Ermakov system for two-dimensional time-dependent 
Bose-Einstein condensates do indeed possess Noether point symmetry. The 
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associated Noether constant of motion and the Ermakov invariant can then 
be used to construct the exact solution for the equations of motion. 

We remark that Bose-Einstein condensation cannot occur in uniform two- 
dimensional and one-dimensional gases at finite temperature, due to destabi- 
lizing thermal effects. However, we are not faced with this problem here since 
we treat non-uniform condensates under the influence of harmonic traps. 
Also, we stress that completely integrable Lagrangian Ermakov systems of 
similar character also arises in the study of the propagation of elliptic Gaus- 
sian beams in nonlinear optics |23], j24fl . These problems in nonlinear optics 
are also described by a nonlinear Schodinger equation formally identical to 
the Gross-Pitaevskii equation. However, the distinctive feature of the present 
work is the presence of explicitly time-dependent harmonic traps, as far as 
we know a feature not previously considered in the literature p3[ p4fl . As 



will be seen in Section 4, there are concrete examples in which the details of 
the time-dependence of the harmonic fields can decisively influence on such 
relevant questions such as the critical time for collapse of the condensate, 
when collapse is in order. 

The work is organized as follows. In Section 2, we describe the basic steps 
for the conversion of the problem of solving the Gross-Pitaevskii equation in 
two spatial dimensions to the problem of solving a coupled, nonlinear system 
of second-order ordinary differential equations. This reduction can occur 
in view of the Gaussian ansatz for the condensate wave function. Thanks 
to the harmonic character of the trap, we formulate the reduced system 
as an Lagrangian Ermakov system. In Section 3, we use the Lagrangian 
formalism to derive an additional invariant by means of Noether's theorem. 
In Section 4 we present the main physical consequences which follow from this 
methodology. In Section 5, the joint conservation of the Ermakov invariant 
and the Noether invariant allows to obtain the exact solution of the system 
in terms of a particular solution for an auxiliary equation whose form is to be 
considered in the following. The whole formalism is exemplified in Section 
6 for some particular choices of time- dependent frequencies of the harmonic 
trap. In Section 7 we present our final remarks and conclusion. 



3 



2 Bose-Einstein Condensates as Ermakov Sys- 
tems 



Our starting point is the Gross-Pitaevskii equation ]T|, 

-^vV + <#lV + £^ = ^, (l) 

2m ot 

which gives a mean-field description for Bose-Einstein condensates. Here, 
ip = -?/>(x, t) is the condensate wave function, while g is the coupling constant. 
In addition, m is the mass of the particles of the condensate and U = t/(x, t) 
is a confining potential. More specifically, we are looking for the two-dimen- 
sional case so that ip = ip(x, y, t), U = U(x,y,t). 

In what follows we will consider time-dependent harmonic traps. More 
specifically, we put 

U = ^l(t){x> + y 2 ) . (2) 

The two-dimensional Gross-Pitaevskii equation can be derived from 
the variational principle 5 S = 0, with 

s = fw{i% + ^-u-iw)4, (3) 

where <ix = dxdy. 

Following references ||-f|, we introduce a Gaussian wave function, 

/ N \V2 x 2 2 

i> = ex p(-^ - e MKPix 2 + foy 2 )) , (4) 

where a, and fa are real functions, depending only on time, and iV is the 
total number of particles. The trial wave function is normalized as 

jW = N. (5) 

The Gaussian ansatz is a reasonable proposition for weak coupling since the 
ground state of the simple harmonic oscillator is a Gaussian wave function. 
Notice, however, that we are admitting a time-dependent frequency u±(t) 
here. 
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After inserting the Gaussian ansatz on the two-dimensional Gross-Pi- 
taevskii action, the result is 



S = J dt£(a,P,j3) , (6) 



for 



N ■ ■ 97) 2 



h 2 1 1 rnw^fy 2 2 gjV 
2m y a\ o? 2 2 v 1 iy 27ra 1 a 2 



Regarding a and /3 as coordinates on configuration space, we obtain the 
following Euler-Lagrange equations for the Lagrangian C, 



and 



ma x ma 2 



2 y 

«2+^ ± (t)«2 = Z^73+T7— — 2— • ( 10 ) 



ft 2 | gN 

m 2 a 2 2irmoi 2 a>i 
Let us introduce some scalings. Let 



Nmujlg . . 

U_l = UJqUJ , t = U t, 9 = - 2 (11) 

2,11 fl 

di = ai/a , a = (— — ) 1/2 . (12) 

mujQ 

Here, Uq is a typical frequency, for instance the temporal average of ui(t), 
and ao is the oscillator typical length. Also, the quantity g can be shown 
to correspond to the ratio of the mean potential energy to the mean kinetic 
energy. 

Dropping bars, the result is 

1 



ai+uj (t)oti = - — 2 F(a 2 /a 1 ) , (13) 



■' 1 1 }{\\ = 

a 2 a\ 

a 2 + u 2 (t)a 2 = —^—^G(a 1 /a 2 ), (14) 
a±a 2 
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where the functions F and G are defined by 

F(a 2 /a 1 ) = a 2 /an+g, (15) 
G(ai/a 2 ) = ai/a 2 + g. (16) 

System (pT3|-p^) is in the class of Ermakov systems. Moreover, ( |l3| - |i"3]) is 
an example of the still more specialized class of Lagrangian Ermakov sys- 
tems. These two properties, of being an Ermakov system and of admitting 
a Lagrangian formulation, are decisive to understand the dynamics of the 
system, as will be seen in the forthcoming sections. Once the solution for the 
Ermakov system is found, the phase functions j3\ and j3 2 are derived from 
(H), thus completely describing the time-evolution of the Gaussian conden- 
sate wave function. 



3 Noether Point Symmetries and Invariants 



All dynamical systems of the Ermakov form flT3| - |HD have the constant of 
motion 

I = -(aid 2 — a 2 ai) + J F(X)dX + J G(X)d\, (17) 

the Ermakov invariant, whatever the functions F and G. It can be easily 
proven that dl/dt = along trajectories. For our purposes, it is better to 
use K — I + 1 which for the choice (|l5]-0) is given by 

k = l -{r*ey + v{e), (18) 

V(6) = -^- e (l+gsin26). (19) 

where r = (a\ + a^) 1 ' 2 , tan# = a 2 /ai are polar coordinates, which are the 
most convenient for the analysis here. 

A central question for any dynamical system concerns the existence of 
a sufficient number of constants of motion for integrability. Hence, it is 
important to find additional invariants independent of K. To address this 
question, we first remark that equations (|13| - p^ ) with F and G given by 
( |T5] - |T6|) are the Euler-Lagrange equations for the Lagrangian 



L = \{r* + r^)-\u\t)r 2 -^. (20) 
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The Lagrangian structure suggest the use of variational techniques for the 
search of extra constants of motion. Among these techniques, Noether's 
theorem [57] plays a distinctive role due to the physical appealing of the 



method. 

In our work ||, we make an extensive analysis of the Noether point 
symmetries for Lagrangian Ermakov systems. Referring to || for the details, 
we derive the following infinitesimal Noether point symmetry transformation, 

f = r + eppr, 6 = 9, t = t + ep 2 , (21) 

where e is an infinitesimal parameter and p is any particular solution to 
Pinney's [[28[] equation 

p + u 2 {t)p=^. (22) 

In (^2|), c is an arbitrary numerical constant. Taking a positive c prevents 
p from changing sign, which is a convenient property for our purposes. It 
suffices to put c = 1, a convention adopted in what follows. Also, just one 
particular solution for Pinney's equation is sufficient for the application of the 
theory. In the following, we always consider the particular solution selected 
by the initial conditions p(0) = 1, p(0) = 0. 

Again by the results of M, it can be shown that the Noether symmetry 



(|2l| ) is associated to the Noether invariant 



J= 1 -(pr-pr) 2 + ^ + K^. (23) 



It can be readily checked that dJ/dt = along the trajectories of ([13| - |T4]) . 

In reference || it is also shown that the Ermakov invariant is associated 
to a dynamical Noether symmetry, 

f = r + err, 9 = 9 + e(-r 2 9 + t9) , t = t + er , (24) 

where r = r(r, 9, r, 9, t) is an arbitrary function. It is apparent that there is 
no choice of r for which ([24]) can be put in the form of a point symmetry, 
independent of velocities. 

The Noether invariant allows to quickly find the exact solution in terms 
of a particular solution of Pinney's equation. Let the quasi-invariance trans- 
formation 

R = r/p, T = [ d\/p(X) 2 . (25) 
Jo 
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In the new coordinates (R, T), the Noether invariant becomes the energy-like 
quantity 

where 

R 2 K 

w{R) = y + W (27) 

We now proceed to the main physical consequences. Integration of the system 
Is given in Section 5. 



4 Physical Interpretation 

The Noether invariant is formally the energy function of a one- dimensional 
time-independent singular oscillator, whose trajectories can be found by 
quadrature. Referring to the potential V(R), we know in advance that 

K < (28) 

is a necessary and sufficient condition for R(T) to attain the zero value at 
some critical rescaled time T = T*. Indeed, for K > 0, there is a repulsive 
force preventing collapse. For K = 0, V(R) is the potential for the simple 
harmonic oscillator, which goes to the origin at some time. Finally, for K < 0, 
there is an attractive force which constrain R to be zero at some time. 

The main physical consequence of our approach is the following. Exam- 
ining the quasi-invariance transformation fl25|) and taking into account that 
p{t) does not change sign (due to the properties of Pinney's equation), we 
immediately conclude that ( p8"D is also a necessary and sufficient condition 
for r(t) to attain the zero value at some critical time t = t* . At this criti- 
cal time, the density of the gas becomes infinity, characterizing the collapse 
of the condensate. At t — t*, the Gross-Pitaevskii equation ceases to be a 
reasonable model, since too great densities are clearly not acceptable in the 
framework of mean field theories. The form of t* depends on the details of 
the frequency u>(t), but in all cases the ultimate evolution of the system is 
dictated by the value of K, which is essentially the Ermakov invariant for 
the problem. 

Looking more carefully the invariant K as given by (|TS|), we see that 
1 + g sin 20(0) < is a necessary condition for K < 0. As «i(0) and «2(0) 
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are positive by definition, we have sin 26(0) > 0. Therefore, we conclude that 

9 £ -sdsjoj (29) 

is a necessary but not sufficient condition for collapse. Indeed, K is composed 
by a positive definite part, related to the angular momentum of the conden- 
sate, and a term which can be non-positive when (^) is satisfied. Therefore, 
attracting Bose-Einstein condensates with negative g satisfying (p9|) does not 
necessarily collapse if the initial angular momentum is sufficiently high. This 
is a distinctive feature from previous works in which the isotropy condi- 
tion a\ = a% sets the angular momentum to zero. For the isotropic cases, 
we also have sin 29 = 1 and K = 2(1 + g), so that the condition for collapse 
is simply g < — 1. To the best of our knowledge, the fundamental role of 
angular momentum in the collapsing or not of two-dimensional Bose-Einstein 
condensates does not seem to be previously pointed out in the literature. 

We now provide some details of the mathematical analysis. Our main 
result is equation (^) bellow which gives the critical time t*. 



5 Integration of the System 

As for any one-dimensional time-independent potential system, the trajecto- 
ries in rescaled variables can be obtained from the second integration of the 
energy-like form (p6|), 

1 r R d\ 

T + ki, (30) 



V2J (J-W(X)) 1 / 2 

where ki is a numerical constant. For simplicity only, we set arbitrarily the 
initial condition as -R(O) = 1, R'{0) = 0. For this choice fl3"0| ) implies 

R 2 = cos 2 T + 2K sin 2 T , (31) 

showing in a more explicit way the role of the sign of K in the existence or 
not of collapse. Indeed, for K > 0, R 2 is a positive definite quantity for all 
T. 

For non-positive K, the rescaled critical time for collapse is given, from 
R(T*) = and (0), by 



T* = arctan . (32) 
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The physical critical time t* follows from 

-■** dX 



o p 2 (A) 




arctan —== . (33) 



The physical critical time strongly depends on the form of p(t) solving Pin- 
ney's equation, and then on the time-dependence of the frequency of the 
trap. 

Expressing (|3~ID in terms of physical coordinates, there results 



r2 -H^^)' +2A - sm2( i'^)') p2(() - (34) 

which furnishes the exact solution for the radial variable in terms of a par- 
ticular solution for the Pinney equation (|22|) . This particular solution may 
be obtained numerically whenever necessary. The exact solution ( |34"D is such 
that r(0) = 1, f (0) = 0. 

So far only the time evolution of the radial variable was considered. In 
order to obtain the angular variable, we integrate a second time the Ermakov 
invariant as given by (0), 

1 r e dX ^ dX 

+ k 2 , (35) 



J {K-V(X)) 1 / 2 J r 2 (A) 



where k% is a numerical constant. The evaluation of the right-hand side of (35) 



depends on the form of r(i) and then on the form of p(t). The left-hand side 



can be evaluated in terms of elliptic integrals |23[ . In the continuation, we are 
mainly concerned with the time-evolution of the radial variable. However, 
notice that a consequence of (^) is that collapse, when it occurs, take place 
with angular velocity going to infinity |23| . 

In the next section, we illustrate the theory for some specific forms of the 
frequency. 

6 Examples 

6.1 Constant Frequency 

In the case of constant frequency, we can set 

uj = 1 (36) 
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in our non-dimensional variables. Then, the solution for Pinney's equation 
with p(0) = 1, p(0) = is 

p = l. (37) 

Hence, the quasi-invariance transformation (p5|) becomes the identity and 
the formulae fl3T| - |3"2"D translate directly to physical variables. In other words, 

r 2 = cos 2 1 + 2K sin 2 1 , (38) 

and, for non-positive K, collapse take place at 

t* = arctan ( I . (39) 




This critical time is a decaying function of \K\, as it should be. That is, the 
more negative is K, the faster is the collapse. 

6.2 Algebraic Decay of the Frequency 

The frequency 

" = (T^y (40 » 

illustrate the relevant case of a trapping potential that decays at an algebraic 
rate. The numerical coefficients in (|40| ) are chosen so as to obtain simple 
expressions in the following. For (fTO|), the solution for Pinney's equation 
with the prescribed initial condition is 

tt'-e + 1\ 1/2 , v 

The final expression for function of t which follows from ( j34[) is given 

by 

^(l-Q' + Mtf 



showing that r — > t asymptotically. Using (p3|), the critical time for collapse 
when K < is obtained as 



r = ^1 + 1^1/2-^1^1/2, (43) 

again a decaying function of \K\. Other kinds of decaying traps are also 
amenable to exact calculations, as for the case of an exponential decay, which 
can be treated in terms of Bessel functions. 
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6.3 Sudden Jump in the Frequency 

As a final illustrative example, consider the case of a sudden jump in the 
frequency at some time r > 0, 

(0, t <t; 
U= {l, t>r. ^ 

For this choice, the solution for Pinney's equation with the prescribed initial 
condition is 

p f(l+* 2 ) 1/2 , t<r; 

P \(c0s(t-T) 2 + (TC0s(£-T)+sin(t-T)) 2 ) 1/2 , t>T. K ' 

This solution is continuous and has continuous first derivative at t = r. Using 
(|34T), we obtain 



{l + 2Kt 2 ) 1/2 , t<r; 

(cos(t - r) 2 + 2K (r cos(t - r) + sin(t - r)) 2 ) 1 / 2 , t > r 



(46) 



for the time-evolution of r, which becomes to be oscillatory at t = r if it has 
not collapsed before. Indeed, by (|33D the critical time for non-positive K can 
attain two expressions according to the values of r and \K\, 



r + arctan(l/ v /2|ir| - r) , T<l/y2|iT 

No matter the values of the parameters r and \K\, the critical time is a 
decaying function of \K\ as it should be. 



7 Conclusion 

The main conclusion of this work is the inequality (p8|), which is a necessary 
and sufficient condition for collapse. This condition shows that the decisive 
quantity for collapse is the invariant K, the Ermakov invariant apart from an 
additive constant, and not the coupling constant g. Sufficient angular mo- 
ment can prevent collapse, even in the case of highly attractive condensates. 
This conclusion was guided by the description of the Bose-Einstein conden- 
sate in terms of a completely integrable Ermakov system. We remark that 
the exact solution shown in Section 3 is the exact solution for a variational 
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approach based on a Gaussian ansatz. This gives a true exact solution for 
the Gross-Pitaevskii equation only in the Thomas- Fermi limit, when the first 
terms at the right-hand sides of equations (|9f]T0|) are disregarded. Neverthe- 
less, the exact solution is an improvement over isotropic solutions (with zero 
angular moment) or self-similar solutions [p9fl . Furthermore, our techniques 
allow possible time-dependence of the trap. 

A natural question is about the extension of this work to three-dimen- 
sional Bose-Einstein condensates. In this case, the Gaussian ansatz for the 
Gross-Pitaevskii action does not lead to Ermakov systems, but to Ermakov 
systems perturbed by terms directly proportional to the coupling constant g. 
This suggests the use of perturbation techniques for small g. This direction 
is now being pursued in a work in progress. 
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